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Abstract 

A system of one-dimensional Brownian motions (BMs) conditioned never to collide with each 
other is realized as (i) Dyson's BM model, which is a process of eigenvalues of hermitian matrix- 
valued diffusion process in the Gaussian unitary ensemble (GUE), and as (ii) the ^.-transform of 
absorbing BM in a Weyl chamber, where the harmonic function h is the product of differences of 
variables (the Vandermonde determinant). The Karlin- McGregor formula gives determinantal ex- 
pression to the transition probability density of absorbing BM. We show from the Karlin-McGregor 
formula, if the initial state is in the eigenvalue distribution of GUE, the noncolliding BM is a de- 
terminantal process, in the sense that any multitime correlation function is given by a determinant 
specified by a matrix- kernel. By taking appropriate scaling limits, spatially homogeneous and inho- 
mogeneous infinite determinantal processes are derived. We note that the determinantal processes 
related with noncolliding particle systems have a feature in common such that the matrix-kernels 
are expressed using spectral projections of appropriate effective Hamiltonians. On the common 
structure of matrix- kernels, continuity of processes in time is proved and general property of the 
determinantal processes is discussed. 

KEY WORDS: Noncolliding Brownian motion; determinantal processes; random matrix theory; 
Karlin-McGregor formula; multitime correlation fimctions; matrix- kernels; spectral projections. 

1 INTRODUCTION 

In the present paper, we discuss noncolliding Brownian motion (BM) in one-dimension with finite 
number of particles and its infinite particle limits — > oo. The condition imposed to particles 
in the model, not to collide with each other, causes "entropy forces" between all pairs of particles, 
which are repulsive long-ranged interactions proportional to the inverse of distances between particles. 
When we draw sample paths of N particles of the system on the spatio-temporal plane, random sets 
of nonintersecting paths are obtained. Viewing them as random patterns of polymers or phase 
boundaries on a plane, the present system has been used as a model of polymer networks [261 133j . or 
a model showing wetting (or melting) transitions [36] in statistical physics; see also \50\ I92[. Recently 
many authors have reported that notion of noncolliding BM and its discrete counterpart called vicious 
walk J36\ is very useful to analyze the polynuclear growth models [78l HH [82l time-dependent 
correlations of quantum spin chains [TH| , traffic problems |7J , and the Chern-Simons theory [27t [B] • 
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1.1 Observations of Three-Dimensional Bessel Process 



In order to demonstrate the important connection between the random matrix theory [66] and the 
noncolUding BM here we show a couple of observations of the three-dimensional Bessel process. The 
noncolliding BM can be regarded as a multivariate generalization of the three-dimensional Bessel 
process given below. 

Let Bi(t), B2{t), 3^(1) be one-dimensional standard BMs (see Section 2.1 for definition). They are 
assumed to be independent and we consider a 2 x 2 traceless hermitian matrix 

where i = \/— T. Since the four entries (M-^'*(t))i<j.fc<2 are BMs, M(i)(i),i G [0, oo) is regarded as 



a matrix-valued process, which describe a diffusion process in the space of 2 x 2 traceless hermitian 
matrices, which is identified with the three-dimensional real space (R denotes the set of all real 
numbers). At each time t E [0, oo), it will be diagonalized by an appropriate unitary matrix and the 
eigenvalue is given by ibX(t) with 



X{t) = ^{B,{t)Y + {B2{t)Y + {B,{t)Y. (1.2) 

If we consider a Brownian particle in R^, B{t) = {Bi(t), B2{t), B3{t)),t £ [0, oo), the distance of the 
particle from the origin (i.e., the radial coordinate of B{t) ) is given by (II. 2p and thus it equals the 
eigenvalue process X{t) associated with the matrix- valued process M^^\t). This is called the three- 
dimensional Bessel process in probability theory (see, e.g. [ISIITH]), and a simple application of the 
Ito formula gives its stochastic differential equation (SDE) as 

dX{t) = dB{t) + t £ [0, oo), X{0) = x>0, (1.3) 

X{t) 

where B{t) is another one-dimensional standard BM than the above Bj{t),j = 1,2,3. Corresponding 
to the SDE (jl.3p . the backward Kolmogorov (Fokker-Planck) equation for the transition probability 
density pxit, y\x), starting from x > at time t = and arriving at y > at time i > 0, is given by 

d 1 92 Id 

■^px{t,y\x) = 2-Q^Vx{t,y\x) + --^px{t,y\x), (1.4) 
and its solution with limt^Qpx{t,y\x) = 6{x — y) is obtained as 

Px{t,y\x) =^^^Pi,hs{t,y\x) (1.5) 

with 

Pabs(t, y\x) = p{t, y\x) - p{t, y\ - x), 

where h{x) = x and p{t,y\x) is the heat kernel given by ()2.ip in Section 2.1. The reflection principle 
of BM can be used to prove that Pabs(^) y\x) is the transition probability density from x > to y > 
during time t of the absorbing BM, in which an absorbing wall is set at the origin and any particle is 
absorbed if it arrives at the wall. It is a matter of course that h{x) = x is harmonic, d'^h{x) / dx"^ = 0, 
but we dare to say that px is the harmonic transform (/i-transform) of pabs looking at (II. 5|) . We 
introduce another 2x2 matrix 

^ ^' ^ \p{t,y2\xi) Pit,y2\x2)j ^ ' 
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for X = {xi,X2),y = (yi,y2) £ R-^i and consider its determinant 



f2it,y\x) = det \M^^\t,y\x)\. (1.7) 



Then it is easy to see that 



hrt 

Pabs(t, y\oc) = y Y/2(t/2, {-y/2, y/2}\{-x/2, x/2}), x,y > 0. (1.8) 

In summary the three-dimensional Bessel process X{t),t £ [0, oo) has two different reahzations; (i) 
the eigenvalue-process of 2 x 2 hermitian-matrix valued process (jl.ip . and (ii) the /i-transform of the 
absorbing BM with a wall at a; = 0. We also observed that the transition probability density of the 
absorbing BM has a determinantal expression of a 2 x 2 matrix ()1.6p - ()1.8p . If we consider the two- 
dimensional BM, it is represented by motion of a point in the two-dimensional space {xi,X2) G R^. 
We put an absorbing boundary on a line X2 = xi and trace the motion of the point in the region 
W2 = {x = {xi,X2) G : xi < X2}. The transition probability density from x = {xi,X2) £ W2 
to y = (?/i, 2/2) G W2 is generally given by the determinant (II. 7|) . As a special case of it (with a 
time-change t t/2), (jl.Sp is given. 

1.2 Dyson's BM Model, Karlin-McGregor Formula and Noncolliding BM 

The noncolliding BM, X{t) = {Xi{t), X2{t), ■ • • , X]\f{t)),t £ [0,oo) is a conditional diffusion process. 
It has the following two kinds of realizations. 

(i) In order to generate the random matrix ensembles, Dyson introduced NxN matrix- valued diffu- 
sion processes |3U| . For the Gaussian unitary ensemble (GUE), N"^ independent one-dimensional 
BMs are used to assign entries of matrix to satisfy the condition that the matrix is hermitian at 
any time. Eigenvalues are real and define an A^-particle system in one dimension called Dyson's 
BM model (with the parameter j3 = 2 corresponding to GUE). This process solves the SDE (see 
|55j for the proof using generalized Bru's theorem) 

dX,{t) = dBj{t)+ x (t)'^-X (tV l^-?'^^' *e[0,oo), (1.9) 

where B{t) = {Bi{t), ■ ■ ■ ,Bfyf{t)) is an A^-dimensional BM; Bj{t), 1 < j < N, are independent 
one-dimensional standard BMs. The SDE (II. 9p is an A^- variable generalization of the SDE for 
the three-dimensional Bessel process (jl.Sp . It was proved that with probability one Dyson's BM 
is non-colliding |81j . 



(ii) We consider the following subset of R-'^, 

Wat = ja; = (xi, • • ■ ,xn) £ R^ : xi < 2:2 < • • • < XArj. (1-10) 

It is called the Weyl chamber of type ^at-i [HI [55]. The absorbing BM is defined by putting 
absorbing walls at all boundaries of the region W^v, whose transition probability density, when 
the process starts from x £ Wjv at time t = and arrives at y £ W^v at time t > 0, is given by 

fN{t,y\x)= det p{t,yj\xk) , x,y£WN. (1.11) 

l<j,k<N L 

This determinantal expression is known as the Karlin-McGregor formula [30]. (Such a deter- 
minantal formula for nonintersecting paths is known as the Lindstron-Gessel-Viennot formula 
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[621 HO] in the enumerative combinatorics, see also [ZellMlEOlESlEIlllSlEaiZSlIMlEOlE^ 
The noncohiding BM is given by the /i-transform of the absorbing BM [41] 



PN{t,y\x) 



hNix) 



fN{t,y\x) 



(1.12) 



where 



■n{x)= (xk-Xj) = 




(1.13) 



is the Vandermonde determinant [64^ \39\ l89] . 

We note that there appear three different kinds of matrices. The matrices representing Dyson's 
matrix- valued process (eq. ()l.ip for the simplest case), matrices in the Karlin-McGregor determinants 
(eqs. (|1.7p and p. lip ) and that in the Vandermonde determinant (|1.13p . Of course, the equivalence 
of Dyson's BM model (the eigenvalue process of the first kind of matrices) with the noncolliding 
BM implies direct connection between the random matrix theory and stochastic processes. In the 
present paper, however, we will show that the Karlin-McGregor formula is much more important. In 
Section 3 we will show that the Vandermonde determinant appears in the Schur function expansion of 
the Karlin-McGregor determinant. With the combination of these two determinants, the orthogonal 
polynomial method is applicable to study the processes. This method has been developed to analyze 
multi-matrix models in the random matrix theory [66] . 

1.3 Matrix-Kernels and Determinantal Processes 

In Section 2, we will explain that the Hermite orthonormal functions are useful to represent BMs. 
Precise descriptions of facts briefly mentioned above will be given in Section 3. Staring from Karlin- 
McGregor's determinantal expression of transition probability density, we will prove in Section 4 that, 
if we specify the initial configuration as the GUE-eigenvalue distribution, the generating function of 
multitime correlation functions is given by a Fredholm determinant for the noncolliding BM, and 
thus multitime correlation functions are generally given by determinants (Theorem 14. ip . The system 
whose spatial correlations are given by determinants is usually called a determinantal point field (or a 
determinantal point process) in probability theory [85^ I84|. Theorem 14.11 states that the noncolliding 
BM is not only a determinantal point field at any fixed time < t < oo, but it is also a determinantal 
point field on the spatio-temporal plane. We say that the noncolliding BM is a (finite) determinantal 
process to express this situation. 



Determinantal processes are generally determined by their matrix- kernels (see, for example, [96] ) . 

The matrix-kernel of our finite noncolliding BM is expressed by the Hermite orthonormal functions 
{^Pk}'kLo? which is called the extended Hermite kernel in [96] . In Section 5 we will show that the 
asymptotic properties of {^k}kLo completely determine infinite particle limits of the matrix-kernel 
and then the determinantal process. Appropriate scaling limits are performed and two infinite particle 
systems are derived from the noncolliding BM. One of them is the spatially homogeneous infinite 
determinantal process with matrix-kernel expressed by trigonometric functions (Theorem 15. Ij) and 
another is the spatially inhomogeneous one with matrix-kernel expressed by Airy functions (Theorem 
15. 2p . The former kernel is called the extended sine kernel and the latter the extended Airy kernel in 
[96] . We will claim in Section 6 that these three determinantal processes (one finite and two infinite 
systems) and others reported in references [THl HHJ [95l [96l [SH [57] have a common structure; the 
matrix-kernels are expressed by spectral projections associated with appropriate self-adjoint operators 
(effective Hamiltonians) |78] . As explained by Spohn [88] and by Prahofer and Spohn [78j . this common 
feature is shared also with the 1-1-1 dimensional Fermi field in quantum mechanics (see also [35\ 143]). 
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It may be due to the similarity between the Karhn-McGregor formula for noncohiding systems and 
the Slater determinant for free fermion systems with the Fermi exclusion principle [TT]. For finite 
and infinite determinantal processes with matrix-kernels associated with spectral projections, we will 
prove that the determinantal processes are continuous in time (Lemma 17. ip and discuss the bilinear 
forms derived from correlation functions (Proposition I7.2( ). Future problems are given in Section 8. 



2 BROWNIAN MOTION AND HERMITE POLYNOMIALS 



2.1 Diffusion Equation and Hamiltonian of Harmonic Oscillator 

Let {i},J-',P) be the probability space. One-dimensional standard BM starting from a point G R 
is defined as a real-valued stochastic process {B{t,uj) : t G [0,cxd)}, which satisfies the following 
conditions (see, for example, [481 142j). Here a; is a label on sample path, to £ ^l. 

1. B{0,uj) = xq with probability 1. 

2. For any fixed to € 0,, B[t) is a continuous real- function of t with probability 1. (With this 
property we say that the paths are continuous in time.) 

3. For any series of times to = < ti < • • • < ^m, M = 1, 2, • • •, {B{tm+i) - -B(tm)}m=o,i,- -,M-i 
are independent and they are normally distributed with mean and variance tm+i — tm- 

Then if we introduce an integral kernel (Gaussian kernel) 



the probability that the BM stays in an interval [am,^m,]; — oo < am < bm < oo, at each time t. 
m = 1, 2, • • • , M, is given by 



That is, the transition probability density of the BM is given by (j2.ip . Since it satisfies the one- 
dimensional diffusion equation (heat equation) 




(2.1) 






(2.2) 



with the initial condition u{0,x) = 6{x — x'), it is specially called the heat kernel. 
We consider the following transformation of variables, (t, x) i— > (r, C) 



r = r(t) = logt, C = C{t,x) 



X 



(2.3) 



V2t 



and set 



u{t,x) =e^^^+'^'^/^U{T,C)■ 



{2A) 



Then the diffusion equation l\2.2h is transformed to 




(2.5) 
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where 



2de ' 2^ • ^'-^^ 

Note that TCu is identified with the Hamiltonian in the coordinate(C)-i'epresentation of the one- 
dimensional harmonic oscillator in quantum mechanics, if we set the mass of the oscillator m = 1, the 
circular frequency uj = 1, and h = 1. 

Following Dirac's description |28t we consider the real- valued Hilbert space with basis {\Q : 
C G R}, which is orthonormal (C|C') = ^{C ~ C) ^^^^ complete 



/ dC\C){C\ = i. 



(2.7) 



Let Tiu be the operator such that (C'I^hIC) = '^(C' ~ O^H with (j2.6p . We consider a state vector 
l^'(r)), which follows the equation 



d_ 

dr 



^ir)) = -Unu-\]\^{r)). 



(2. 



If we multiply {(\ to ()2.8p from the left, we will have the equation p.Sp with U (r, () = Given 
|^(r')),r' E (—00,00), the solution of ()2.8p for r > r' is obtained as 

Mr)) = exp |-i (hh O} \^{r')). 

Then 

U{t,0 = (C|exp|-i(wH-^)(r-r')}|^(r')) 

^ dC (CI exp |-i (t^h - ^) - \C)U{t', a 

where ()2.7p was used. Insert this result into (12. 4|) . we have 

u{t,x) = I dC e 
J —00 

/oo 
(ia;'p(t - t',a;|x')n(t',3;'), 
-00 

where 



'.-(^K^)/2+((.')+(C')^)/2^^| |-i (un - i') _ r') \ \C')u{t\x' 



1 



-2:'*/4t+(x')"'/4i' 



expi f^H- ^ ) (r-r) 



/2i \ V2t 

This is the transition probability density previously given as ()2.ip . 



/2t' 



(2.9) 
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2.2 Hermite Polynomials and Equalities 

Let N = {1, 2, 3, • • •} and Nq = N U {0}. The eigenvalues of are n + 1/2 with n £ Nq; Wnln) = 
(n + 1/2) Here {|n) : n € Nq} denotes the set of eigenvectors of Tin, which is orthonormal 
{n\n') = dn,n' and complete Yl'^=o = 1- Let {Hn{x) : n G Nq} be the Hermite polynomials 

2 / a \ ^2 



n! 



where denotes the largest number not greater than z. They are orthogonal 

/ dxe~'^^Hn{x)Hm{x) = hn5n,m, n, m G Nq 
Jr 

with hn = A/vr2"'n!. We set 



MO = -^e-^'/2F„(C). (2.11) 



Then {<^n{C) ■ ^ Nq} are orthonormal and 

(C|n) = (n|C) = MC) 



is established. 
Now we define 

Then 



= 7^ ( - 7^ ) . (2.12) 



n 
2 

and (|2.9p gives 



n^\n) = -\n), nGNo (2.13) 



n=0 

^ oo 
= ^e-^/2-CV2+{C')V2 ^^^|g-nr/2|^^^^|gnr72|^/^ 

^ oo 

= i=e-/2-CV2+(C')V2 ^ e-(--')/Vn(C)9^n(C'), 



n=0 

that is 



n=0 

This expression for the heat kernel (j2.1|) is called Mehler's formula [10^ I91j. 
We introduce the vectors 

\t,x) = e^'l^e^^^\C), 



= ^e-(-+f')/2(^|e-H^r_ (2.15) 
v2 
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It is easy to see that orthonornality and completeness of {\Q : C G R-} imply {t,x\t,x') = 6{x — x') 
and 



dx\t,x){t,x\ = l. (2.16) 

R 

By (j2.16p . we have the equalities for < ti < t2 <J ^3 

/ dX2{h,X:i\t2,X2){t2,X2\tl,Xi) = {t3,X:i\ti,Xi), Xi,X3GR, 

/ dxi{t2,X2\ti,xi){ti,xi\n) = {t2,X2\n), n e Nq, X2 € R, 
JR 

dxi / dx2{n\t2,X2){t2,X2\ti,xi){ti,xi\n) = {n\n) = (5.„_„/, n,n' G Nq. 
R JR 



Though these seem to be trivial, if we note that ()2.9p is rewritten as 

p{t-t' ,x\x') = {t,x\t' ,x'), < t' < i, G R, (2.17) 
they become meaningful; for < ti < t2 <^ ^3 

/oo 
(ix2j'(t3 - i2,a;3k2)p(i2 - ii,a;2|a:^i) = p(t3 - *i,a;3|a^i), xi,j;3GR, (2.18) 
-oo 

/oo 
(ixip(t2 - ti,a;2|xi)(/)„(ti,xi) = (/)n(t2,a::2), nGNo,X2GR, (2.19) 
-oo 

/oo /'OO 
dxi / dx2^„(t2,a;2)p(*2 a;2|a;i)(/>n'(*i) 2^1) = '^n.,n', e No, (2.20) 
-00 J —00 



where 



Ut,x) ^ (t,x|n) = i=t-("+i)/V-^/v(^), (2.21) 



^n{t,x) ^ (n|t,x)=r/V'/^Vn(' 



n G No. (2.22) 



The equation (12.180 is called the Chapman-Kolmogorov equation. The equalities p.lOp mean invari- 
ance of the functions (pn^n G No, with respect to the heat kernel. 

3 NONCOLLIDING SYSTMES 

3.1 Application of Karlin-McGregor Formula 

Here we introduce the noncolliding BM in a finite time-period (0, T) with T > 0. It is defined as 
an A^-particle system of one-dimensional standard BMs conditioned not to collide with each other in 
(0, T). As mentioned in Section 1.2, the transition probability density of the absorbing BM in the 
Weyl chamber W^v of type Ajv-i is given by (jl.lip as an application of Karlin-McGregor formula 
[IS]. The probability that B{t) starting from a;' G W^r stays in Wat up to at least time i > is given 
by 

MN{t,x')= fN{t,x\x')dx, x' eWN. (3.1) 
Jwn 

The transition probability density of noncolliding BM is then given by 

ffAr,T(t', x'; t, X) = ^r'''^^~l;% fN{t - t', x\x') (3.2) 
J\In[-L — t , x'j 
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for < i' < t < T, x^x' G Wat [52l [531 El]- It should be noted that this process is in general 
temporally inhomogeneous. In the following, we will consider the T — > oo limit to make the process 
be homogeneous in time. 



3.2 Schur Function Expansion 

By multilinearity of determinant p. lip with ()2.ip . 



j]^(t,x\x') 



1 \ 

— ) e-(l^l'+l^'l')/2* det {e-^^lX 



Yl'j=i^'j- Consider Fj\f{x,y) = deti<jjc<7v[e^^^''] for a pair of multivariates x = {xj}jLi 



where |a;p = 

and y = {yj}jLi- By definition of determinant, F{x,y) is skew-symmetric under any exchange of 
indices of {xj} and also it is for {uj}; Let Sn be the symmetric group of variables (the set of 
all permutations of N variables) and for a £ Sn write (t{x) = (^^(i), • • • , a;o.(jv))- Then for any 
a G Sn Fi\j{a{x),y) = F]\f{x, a{y)) = sgn{a) Fjy [x , y) . A fundamental skew-symmetric polynomials 
of multivariate x is given by a product of differences (the Vandermonde determinant) (|1.13p . The 
quotient of Fjs[{x,y) divided by hN{x)hfyf{y) is a symmetric function both of x and y. The following 
lemma shows an expansion of the symmetric part using the Schur functions, which are labeled by 
partitions fi = {ni, ^2, • • •), sets of nonnegative integers in decreasing order fii > IJ,2 > • • and defined 
by 



det [xf+^-'^l 

l<j,k<N ^ 



det |x, 

l<j,k<N ■' 



N-ki ' 



(3.3) 



The non-zero /Xj's in a partition n are called parts of n and the number of parts is called length 
of fj, and denoted by £{ii) [Ml [39l [89]. The Schur function expansion is a special case of character 
expansions (see [8l [9l [591 [551 [l3] ) . Let T{a) = C^y'^^^dy for a > (the Gamma function) and 
note r(a + 1) = al if a e Nq. 



Lemma 3.1 



det [e^^'^H 

l<j,k<N 



hN{x)hN{y) ^ 
hN{x)hN{y 



Sf,[x)Sf,[y) 



,:ii,)<N ULl nf^k + N-k + l) 



X {1 + 0{\x\)} in \x\ 0. 



(3.4) 



Proof. By multilinearity of determinant, we have 

{xiykY 



det [e'^jy^] 

l<j,k<N 



det 

l<j,k<N 



.71=0 



N 



= V TT ^;^7 r det 

n={ni,n2,---,njv)GN^ "^=1 

We can see that for any symmetric function f{n) of n = (ni, • • • , n^r) G N^, 



{xjykT' 



(3.5) 



fin) det 

^<j,k<N 



[Xjyk) 



1 

m 



(Te<Sjv 



det 

l<j,k<N 



(Xjyk) 
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and 



V det 

^-^ l<j,fc<iV 



det [x^'^l det [y^"]. 

l<j,fc<iV ^ l<i,m<N ^ 



Since deti<j^fc<Ar[xJ*] = if n^^ = n^j for any pair 1 < fci / A;2 < A^, ()3.5p equals 

det «'=) det (y,"'") 

^ l<i,fc<Af l<^,m<Af 

0<ni<n2<-<njv 11^=1 ^C^' + 1) 

Here we change the variables in summation from {nj} to {/^j} by [ij = rij — N + j,l < j < N. Using 
()3.3p we obtain the first equation of (j3.4p . Since s^(0) = unless /i = = (0, 0, • • • ,0) G Nq', and 
so(0) = 1, the estimation in |a;| ^ is given as shown by the second equation of ()3.4p . | 
By this lemma, we have the estimate 

/^(t,^|^') = ^i-A^V2;,^(^)/,^(^')eH^lV2*^|l + o(^)| in ^-0 (3.6) 
with Cn = (2vr)^/2-|-[^^r(j). The integral formula 

Jrn fj; r(l + 7) 



e 



is found in [66] (eq.(17.6.7) p. 321) as a variation of the Selberg integral whose proof was given in 
|63j . If we set 7 = 1/2, a = l/2t and note that the integral over R'^ can be replaced by the integral 
over Wiv multiplied by A^!, we have 

/ e-\''\'/^'hr,ix)dx = C7^t^(^+i)/^ (3.7) 
with = 2^/2 nJLi r(i/2). Similarly by setting 7 = 1 and a = l/2t, we have 

\^\"/^\hN{x)fdx = CNt^'/\ (3.8) 
Using (13. 7p with (13. 6p . we obtain the asymptotics of Mn, 

AA^(t,^') = §^i''^^'^-'^/V(a;')x|l + o(^)| in ^-0. (3.9) 
The integral (13. 8p will be used shortly. 
3.3 Temporally Homogeneous Limit 

By the above estimate (j3.9p . we can take the T ^ 00 limit in (|3.2p and obtain the transition probability 
density, which is homogeneous in time, i.e., a function of time difference t — t', 

PN{t — t' ,x\x') = lim gj\f^T{t' ,x';t,x) 

T — >oo 

= hN{x)fN{t-t',x\x')--^. (3.10) 
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From now on, we consider the noncolliding BM, which is defined by this transition probabihty density. 
It is a temporally homogeneous process and we denote it by X{t) = {Xi{t), X2{t), ■ ■ ■ ,XN{t)). 



Remark 1. The product of differences (the Vandermonde determinant) hi\}{x) given by ()1.13p is a 
harmonic function in the sense 

^ q2 

V^hN{x) = Q^hNix) = 0, 
j=i j 

which has strictly positive values at interior points of W^r and zero at the boundary. Eq. (IS.lOp is 
considered as a transformation from /tv to pn associated with the harmonic function, which is called 
the /i-transform [29]. That is, the temporally homogeneous noncolliding BM is an /i-transform hy 
of the absorbing BM in the Weyl chamber W^v [41j. It is easy to confirm that pAr(t, • \x) satisfies the 
following backward Kolmogorov equation 

d 1 

— ti(t, x) = -V'^u{t, x) + V logh]\f{x) ■ Vu{t, x), 

which is a multivariate extension of (II. 4p . It implies that the process X{t) is a diffusion process, 
which solves the SDE (jl.9p of Dyson's BM model (with (5 = 2), which describes the time-evolution 
of eigenvalues of hermitian matrix- valued (GUE) diffusion process (3m 155]. This equivalence between 
Dyson's BM model for GUE (the eigenvalue part of the matrix-valued BM) and the present noncollid- 
ing BM (BM conditioned not to collide) is a multi-dimensional extension of the equivalence between 
the three-dimensional Bessel process (the radial coordinate of the three-dimensional BM) and the 
one-dimensional BM conditioned to stay in the positive region B{t) > 0, as announced in Section 1. 
See also [53l [58l |55] . 



Let vq^x) be the probability density at time to > 0. Then the probability density of distribution 
X{t),t > to is given by 

vtix) = / pNit — to,x\x')iyQ{x')dx' 
Jwn 

= hN{x) [ fN{t-to,x\x')^^^^dx'. (3.11) 
By the estimate (13. 6p . we will see that 

if the distribution vq has finite moment. Note that the integral formula (13. 8p guarantees that (I3.12p 
is normalized. It should be noted that the distribution (|3.12p is equal to the eigenvalue distribution 
of hermitian random- matrices in GUE with variance a"^ = t [66] . 

Proposition 3.2 For any initial distribution having finite moment, the asymptote of probability den- 
sity of distribution of the noncolliding BM, X{t), in t ^ oo is expressed by the eigenvalue distribution 
of GUE with variance t. 
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4 DETERMINANTAL PROCESS 
4.1 GUE Initial Distribution 

Consider a sequence of times, < to < *i < " " " < ^M, M = 1, 2, • • • for observations of distribution of 
X(t). Given the initial distribution z/q at time to, multitime probability density is given using (|3.1U|) 
as 



M-l 



PN{to,x^'^^-ti,x^^'^;---;tM,x^^'^) = PN{t 



m=0 



M-l 



m=0 



m+1 i^m 



Now we assume that fo is the GUE-eigenvalue distribution with variance to, 

Cat 

Then ()4.ip becomes the product of /at's multiplied by a factor 



(0)^ 



(4.1) 



(4.2) 



(4.3) 



By multilinearity of determinant and the fact that the coefficient of the highest order term of the 
Hermite polynomial Hn{x) is (2x)"' (see (I2.10p ). the following equality is established, 



det 

l<j,fc<Af 

Then if we define the multivariate functions 



det 

l<j,k<N 

det 

l<j,k<N 



(pj-lito,Xk) 
(t>j-l{tM,Xk) 



(4.4) 



where 4>n and 0„ are given by (j2.2ip and (|2.22|) . respectively, the factor (j4.3p is readily shown to be 
equal to fiQ{x^^^)fj,Mix^^^^). That is, the multitime probability density (14. ID is written as 



M-l 



(0)^ 



m=0 



This expression shows that the probability law is invariant under any permutation of indices of 
particles ; Xj™^ i— > x^^^-^ , < m < M, a G Sj\[ ■ Then description will be easier if we regard that 
particles are identical and indistinguishable. We denote by X the space of countable subsets ^ of R 
satisfying ti(? H K) < oo for any compact subset K. The space X is a Polish space with the vague 
topology. For x = {xi,X2, ■ ■ ■ ,Xn) G U^i^^i denote {xj}^^^ S X simply by {x}. Then we 
consider the process on the set X, H^(t) = {X{t)},t G [to,oo), such that, for any M + 1 times 
to < ti < • • • < Im-i < thi < oo (M = 0, 1, 2, • • •) the multitime probability density is given of the 
form 



Pn [to, {^(°)}; ti, {^^^n; • • • ; tM-i, {x("'-'>};tM, {x^"'>} 
m{{x'^^'^]) n /iv(Wi-t™,a.("+')|^("^Vo({a;(°)}), 

m=0 



,(A^-l)l.- 



,{My 



(4.5) 
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For aj^"^) E R^, < m < M, and iV' = 1, 2, . . . , N, we put tcj^^ = (^j^S'"^ x^"'^ . . . , j;}^^) . For a 

sequence {A'ml^^Q of positive integers less than or equal to A^, we define the (A'^o, . . . , A^M)-iiiultitime 
correlation function by 



P;v(to,{^(°^};...;tM,{^(^')}) n 



M N 

(N-N^)l . ii • 



(4.6) 



Expectations with respect to the configurations {X(to)}, {^(h)}, ■ ■ ■ , {^{tM)} are denoted by Eat : 
Eat 



f{{X{to)},{X{h)},...,{X{tM)}) 

M+l 



1 \ ^"+^ r / A/ TV 

^) / /({.(°)},...,{.w})p^(to,{-(°)};...;iM,{-(^^^}) n n^-r- 

(4.7) 

Remark 2. Set t' = 0,t = to > 0,x = a^^^) in (fXTO]) and consider the limit x' 0. By (fSlB]) . we 

can see 

hm p^(to,a;(o)|a;') = ^^o(^(°^). (4.8) 
|a;'|-^o 

This fact is essentially the same thing with that stated as Proposition 13.21 for the scaling property of 
the process. The GUE-eigenvalue distribution ()4.2p adopted here as the initial distribution at time 
io > is immediately realized if we start the present noncolliding system at time from {0}. The 
state {0}, which is a boundary of W^r, is entrance [52l [SSj [5i] . 

4.2 Generating Function 

Let Co(R) be the set of all continuous real functions with compact supports. For / = (/o, /i, • " " i /a/) G 
Co(R)^^, and = {6q,9i, ■ ■ ■ ,6m) G R*"'^, the generating function for multitime correlation functions 
is defined for the process {X{t)},t £ [0,T] as 



*7v(/;0) = E 



N 



M N 

exp {J2^-^Y1 fm{Xj^{tm)) 
m=0 jm=l 



(4.9) 



Let 



Xm{x) = e^'"-^™'^'') - 1, < m < M, 
and write (j4.9p as ^'Ar[x]- Then by the definition of multitime correlation function ()4.6p . we have 



N N 



(0)i.. r^(i)i. 



R™i 



(1) 



TV M 

E E - E n^^ 

No=ONi=0 NM=Om=0 
M 

^ n Ylxm^x^J^^^pN [to,{x'-^l};ti,{x^^l};...-tM,WNj}) ■ 

m=0 j=l 



Nm 



(M) 



(4.10) 
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By the definition (fOj) with ([iTT]) and (f^ . we have 



^n[x] 



1 

iV! 

M-l 



M+1 



M N 



n n^-r . det 0,-i(tM,xr)(i+xM(xr)) 



RiV(Af+l) 



r?i=0 j=l 



l<j,k<N 



l<j,k<N L 



m=0 



det 

l<j,k<N L 



(/)j_l(t0,4°'') 



By definition of determinant, it is easy to prove the following identity for square integrable continuous 
functions gj,gj,l < j < N, 

N . 

Y\dxj= det / gj{x)gk{x)dx 
which is called the Heine identity. By repeated applications of this identity we have 



T7J I det 

N\ JyiN l<j,k<N 



9j[Xk) 



det 

l<j,k<N 



9j[Xk) 



(4.11) 



(4.12) 



with 



M 



Fjkix] = I ndx("){^,-_i(tM,x(*^))(l + XM(x(*^)))| 



m — i 

J] {p(Wl-^™,x('-+l)|x(™))(l + x„^(x('")))}0fc_l(^o,: 



m=0 

M-l 

X ii {p{t 

m=0 

By the Chapman-Kolmogorov equation (|2.18|) and the invariance (|2.20p . i*'jfc[0] = 6jk and then ^'7v[0] = 
1, which implies that ()4.5p is indeed normalized. If we use the notations introduced in Section 2, it is 
written as 

M 

Hdx^^^ 0-_l|t,,,^W){l+;,^,(xW)} 



Fjkix] 



m=0 



M-l 

X J] [(Wi,x^"^+'^|im,rE(™)){l + Xm(x('"))}](to,x(°)|A;-l) 

= {j-l\k-l) 

M+1 „ e 

+ E E / ^n^^^""^0--l|tn^i,X^"^^^)Xmi(x("^^))(t„.,,x("^i)|t„„:E(-^)) 

i=0 A4>mi> - >mi>0 ^' n=l 

X • • • X {t-iYi^_-^ , X^ * , X^ Xm^ (^^^ ^"^ ) (^m^ 1). 

Now we introduce an indicator 1j_ such that 



{tm,x\tn,y) if m>n 
otherwise. 



(4.13) 
(4.14) 



{tm,x\l+\tn,y) = 

Then ()4.13p can be written of the form of an infinite series 

Fjk[x] = (j-l|A:-l) 

oo M M e 

+ E E ••• E / ^n^^^""^o--iiv,^^"^^^>x™i(x^'"^^)(v,x('"^)|i+it„,,x(™^)) 

£=Omi=0 m^=0 n=l 

X • • • X (tm^_ J , X^ ^ I 1 , X^ ^'^)Xmi ^"^ ) {trrii i X^ ^"^Ik 1), 
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in which all terms with i > M + 1 are zero by ()4.14p . Let 

X+'''{x,y) = {tm,x\i+\tn,y)Xn{y), m,n e {0,1,- ■■ ,M}, x,y (4.15) 

and define an M x M matrix- kernel x+i whose (m, n)-element is given by the integral kernel ()4.15|) . 
The (m, n)-element of its q-th power {x+)'^y = 2, 3, • • • , will be the following kernel given by {q — 1)- 
multiple integral, 

A/ „ A/ „ 



mq-i=0 ' 



Then we have 
FMx] = {j-l\k-l) 

M „ M „ oo 

+ E / ^^E / dy(j-l|t„,x)xr„(x){5„,„5(x-y) + J^[(x+rr'"(^,y)}(in,y|A^-l). 

m=0 n=0 g=l 

(4.16) 

Let l"*'"(x,y) = (^m,n'^(a; — y). Since it is easy to see that 



A/ „ oo 

Y / ciy[i-x+]™''(x,y){5,,n%-^) + 5^[(xH-)t'"(y,^)} 

£=0 •'^ 9=1 



i'"'"(x,z) 



(4.17) 



we can write 



5^,J{x-y) + Y[{x+rr'^{x,y) 

q=l 



1 - x+ 



ix,y). 



Then (I4.16P becomes 



where 



i?^Vx) 



M 

E / ^-^ ~ i|*n,y)xn(y) 



(4.18) 



1-X+ 



(y,a;) 



Cl."'^(^) = (tm,x|fc-l). 



(4.19) 



4.3 Predholm Determinant and Determinantal Process 

Let Bf'\x) and c'^"'\x), < m < M, 1 < j < iV, be square integrable continuous functions. Then 
the following formulae can be proved; 



M 



det 

l<j,k<N 



(m) / 



Y / dxBp{x)Clr\ 



det 

l<j'-"^\j''"^ <N,0<m,n<M 



+ / dxB|^^)(x)c|™}(j;) 



R 



iV iV 



Af M 



Aro=OAfi=0 Arjv^=Om=0 ™ ^ " j=l j=l 



.(1) 



\[dx 



(M) 



N 



X det [Vcf)(x;"))4")(xl")) 

p=l 



(4.20) 
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The expansion formula of the last expression defines the Fredholm determinant, which is abbreviated 
as 

N 

i-'"(rr,y)+^Cf)(x)5(")(y) 
p=i 



Det 



The generating function (j4.12p with (|4.18|) and (j4.19p is then expressed as the Fredholm determinant, 



M „ 

+ ^ / dy{t^,x\V''\te,y)xi{y) 



where 

Note that (filTl) implies that 

M 



N 



(4.21) 
(4.22) 



i-'"(x,z)=f] / dy{i^''{x,y) - {t^,x\U\t,,y)xi{y)} 



Plugging this into (|4.2ip , we have 



^n[x] = Det 



M 



£=0 



/ dy[i^'\x,y) + {trn,x\{r^ -i+)\ti,y)xi{y) 



(y^z) 



Det[i'"'"(x,y) + {tm,x\{V'' - U)\t„,,y)xn{y) 



xDet 



{x,y) 



Det 



l-'-(x,y) + {tm,x\iV'' - l+)\tn,y)Xniy) 



(4.23) 



Here 



Det 



1-X+ 



m,n 



ix,y) 



l/Det 



[i-X+r'^{x,y) 



and we have used Det [1— x-|-]'"'"(x, y) = 1, which is concluded from the fact that [1— x+]™''"(x, y) = 

for m < n by the definition of x+; (|4.15p with ()4.14p . and it is 5{x — y) for m = n. 
Let 



S^'^^ix^y) = {tm,x\V \tn,y) 



Sn' ix,y) = -{tm,x\{l+-V )\tn,y)- 

Following the formulae given in Section 2, S^'"" is determined as 

N 

S^'^^ix^y) = Y{tm,x\p-l){p-l\tn,y) 
p=l 
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1 



N 

.g-a;2/4tm+s/2/4tn ^ / ^ 

p=l 

N-1 

' fc=0 



(p-l)/2 



fc/2 



1 ) (p- 1 



Combination of this with (j2.14p gives S^'^ as 

5^7v''^(a;; y) = ^'^'"'{x, y) - l{m>n}P{im - tn,x\y) 



1 



N-l 



1 



fc/2 



1 n \^rn 



X 



fc=Af 



Vk 



Vk 



if m <n 
if m > n, 
otherwise. 



(4.24) 



where Ij^} is the indicator of a condition uj; = 1 if is satisfied and 

As shown above, not {S^'"'{x,y)} themselves, but determinants of matrices made of them are 
observables. By definition of determinant, factors {e~^^/^*'"+2/^/4tn| q£ in (I4.24p are 

completely cancelled out, when we calculate determinants. So here we define the following matrix- 
kernel by omitting these factors in {S^'^{x,y)}, 



KN{tm,x;tn,y) = < 



^-1 / ^ \k/2 



k=0 



^k 



^^"^ fc=0 

oo / ^ s^k/2 



Vk 



k=N 



tr 



^k 



"Pk 



y 



if m <n 
if m > n. 



(4.25) 



and rewrite (|4.23p as 

y)xn{y) 

This Fredholm determinant is by definition expanded as 

Det [^(^m.n'^Ca; -y) + ^N{tm, x; tn, y)xn{y) 



(4.26) 



N N N M „ No ^ A^i 



ATq =0 Afj =0 Nm=0 rn=0 

M N,n 

J™) 



Nm 
Hdx 

i=i 



(M) 



X 



m=0 j=l 



det 

l<j<Nm,l<k<Nn,0<m,n<M 



(4.27) 



Comparison of (|4.27p and (j4.10p determines all of the multitime correlation functions. Now we sum- 
marize the above results theorem. 

Theorem 4.1 The temporally homogeneous noncolliding BM, H^(t) = {X{t)}, starting from the 
GUE- eigenvalue distribution ((^.^[ ) at time to > 0, is a finite determinantal process in the following 
sense. 

(i) The multitime generating function is given by the Fredholm determinant ^J^^E), where the 
matrix-kernel K^r is given by \4.25^ . 
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(ii) Any multitime correlation function is given by a determinant; for any M > 0, any sequence 
{^m}m=o of positive integers less than or equal to N, any time sequence to < ti < ■ ■ ■ < tM < oo, 
the [Nq, . . . , Nm) -multitime correlation function is given by 



det 

l<j<Njn,l<k<Nn,0<m,n<M 



KN{tm,xf^^;tn,X^l^^) 



(4.28) 



The following relations hold for the Hermite polynomials, 

Hk+i{x) = 2xHk{x) - 2kHk-i 



(4.29) 
(4.30) 



From 



^Hk{x) = 2kHk^i{x), A; = 1,2,3, 
(|4.29p . the Christoffel-Darboux formula is derived for the Hermite orthonormal functions {^k{x)}k&'t^o > 



JV-l 



V 2 X — y 

for X ^ y. Eq. ()4.30p can be used to evaluate the limit y ^ x in ()4.3ip and we find 

N-l 2 2 

^[^k{x)] =N^LPn{x)^ - V^i^ + l)^N+l{x)ipN-l{x). 
k=0 



if m = n, 



Then the matrix-kernel have the following simpler expressions 

fW ipN{x/V2t^)ipN^i{y/y/2t^) - ipN~l{x/y/2t^)ipN{y/y/2t^) 

V 2 x-y 

if X ^ y 

1 ( ( X nrmr; — ( x 



'^N{tva-,X\tm-,y) = < 



/2tr 



if X = y. 



(4.31) 



X 



(4.32) 



5 INFINITE PARTICLE SYSTEMS 

5.1 Wigner's Semicircle Law 

The density of H^(t) is given by 

7V-1 



PNit,x) 



as a special case (M = 0, A'^o = 1 with setting to = t) of Theorem 14. 1 1 with (j4.32p . It is easy to confirm 
that PN{t, x)dx = N hy the orthonormality of Lpk{x). The following estimations for asymptote in 
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N ^ oo are established [101 191j . Let e and uj be the fixed positive numbers. We have 

1/4 



V / vvrsmc/) \iV 



X < sm 



^ + ^') (sin20-20) + ^7r 



(m) v^at ( \/ 2 + 1 cosh ( 



1 



1 



V27r sinh 4> \2N 



1/4 



X exp 



Y + (2</.-sinh20)+^7r 



l + O 



N 



Using them, we wiU have the asymptote of the density profile in — > oo, 

1 



PNit,x) ~ < 



'2N- — 
TrV2t V 2t 







e < 6 < UJ. 



if 



2VNt <x < 2VNt 



otherwise. 



(5.1) 



The distribution of N particles has a finite support, whose interval oc v -/V, and thus /9 ~ v — > oo as 
N ^ CO for fixed < t < oo. If we set x = 2-\J Nt^, we see 



lim —PN(t,x)dx= < 



-\/\-i^di if -\<i<\ 

vr 







otherwise, 



(5.2) 



which is known as Wigner's semicircle law [66]. See also jSlJ. In the following we consider scaling 
limit, in which long-term limit t ^ oo is taken at the same time with — > oo. 



5.2 Bulk Scaling Limit and Homogeneous Infinite System 

First we consider the central region x ~ in the semicircle-shaped profile of particle density in the 
scaling limit 

t ~ AT ^ oo. (5.3) 

In this limit the system becomes homogeneous also in space with a constant density p = l/vr. We call 
this the bulk scaling limit. 

Theorem 5.1 For any M > 0, any sequence {Njn}m=o of positive integers, and any strictly increasing 
sequence {sm}m=i of positive numbers 



lim pn{N, {xPJ; N + 2si, ■■■■,N + 2sm, {x^^'J}) 



det 

l<j<Nm,l<k<N„fl<m,n<M 



(5.4) 
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where 



- / (iue('*~*)"%os(u(y - x)) if t<s 
71" Jo 



K{t,y;s,x) 



sm{y — x) 
Ti{y - x) 



if t = s 



1 

— / dne~(*"^)"'cos(n(y-x)) if t 
71" J I 



> s. 



(5.5) 



Proof. For any u E R, the formula 



lim (-l)V/^ 



2^/1 



COS u, 



lim(-l)V/V2£+i 



u 



2Vi 



: smu 



/TT 



are known |1UI [M] . We note that 



~ e 



1 + ^ 

2{Sn~Sm)a/N 



N 



1 + 



a/N 



for ^ 1 with fixed number a. Then (j4.25p with m < n is evaluated in A'^ ^ oo as 



^N{tm,y;tn,x) 

N/2-1 



)e/N I ^ 



COS 



l2i \ 
-X +sm 



- E 

ttN ^ 

1=0 

^ j dAe(''""'*'")'^-^|cos(y\/A)cos(x\/A) +sin(y\/A)sin(x\/A)| 

1 2 

-/ (itie(''""*")"%os(n(y - x)). 
vr ./n 



(5.6) 



\ . / /2£ ^ 



In particular, when m = n, i.e., Sn—Sm = 0, the integration is readily performed to have du cos(u(y- 
x)) = sin(y — x)/{y — x). Similar evaluation in — > oo can be done also for ()4.25p with m > n. § 



In the bulk scaling limit (j5.3p . the temporally and spatially homogeneous infinite particle system 
is obtained, whose multitime correlation functions are given by ()5.4p . The matrix-kernel ()5.5p is called 
the extended sine kernel in [96] . In the present paper, we will call it simply "sine kernel" . The system 
with the sine kernel was studied by Spohn ^88j , Osada |72l [7^ , and Nagao and Forrester [TOj , as an 
infinite particle limit of Dyson's BM model with /3 = 2 [30]. See also [96ll2]. 

5.3 Soft-edge Scaling Limit and Spatially Inhomogeneous Infinite System 

Next we consider the scaling limit 

t~Ar^/3 and x ~ 2A^2/3^ ^^ ^^ 
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Since ()5.7p gives /2t ~ 2A^, the vicinity of the right edge of semicircle-shaped profile ()5.ip will be 
closed up, and we will obtain a spatially inhomogeneous infinite particle system in this scaling limit. 
Following the random matrix theory [66j . we call ()5.7p the soft-edge scaling limit. 
In order to describe the limit, we introduce the Airy function 



Ai(x) 



27r 



dk e 



It is a solution of equation 



which behaves as 



Ai(x) 
Ai(-x) 



-— :jAi(x) = xAi(x) 



2^x1/4 ^\ 2, 



(5.8) 



(5.9) 



vrx 



1/4 



COS I -x^/^ - ^ 



m X ^ oo. 



In the proof of the following theorem, we will use the formula 

lim 2-^4^1/12^^ ("y^ + ^r^A = M{u) for ueB.. 

l^oo \ -^2 / 

Let 

aN{s) = 2iV2/3 ^ 2iVi/3s - s^ 
and a;7V'('S) = {cln{s) + xi, OAr(s) + 2:2, • • • , aAr(s) + xn')- 



(5.10) 
(5.11) 



Theorem 5.2 For ani/ M > 0, any sequence {Nm}^=o of positive integers, and any strictly increasing 
sequence {sm}.m=i of positive numbers 

hm pN{N^lMx%lm-N^'^ + 2s^,{x%\{sr)]- • • • ; ivV3 + 2sm, {x^n^{sm)]) 



det 

l<j<7Vm,l<fc<7V„,0<m,n.<M 



= PAi (0, ^^Q^ ; Si , ^j^^V • • ; SM, ^ 



(M) 



(5.12) 



where 



( /-o 



K:{t,y-s,x) 



d\e^'-^^^ki{y-\)Pd{x-\) if t<s 
dA e-(*-")^Ai(y - A)Ai(x - A) if t> s. 



(5.13) 



Proof. Putting k = N — p — 1 in the summation of (14.250 for m < n, we have 



KN{tm,y;tn,x) = ( 



{N-l)/2 



p=0 



-p/2 



2tn 
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With the same reason mentioned above eq. (I4.25p . we can omit the factor (tn/tm)^^ ^)/'^ , Since, when 
we set tm = N^/^ + 2sm, 



aN{sm) + X 



V2 



we can use the formula (|5.1U|) : 



X 



^ + ^N-^/^x] 
V2 J 

~ ipN-p-i (^2(iV -p-l) + l={N-p- l)-i/6 |: 



P ] 

Ari/3 J 



For 



^ 1 + 2s„/iVV3 1 



-(^fn-S™)p/iVl/3 



in iV 



oo, 



we have 



N-l 



P 



Ai x 



p=0 



)"Ai(y + ■u)Ai(x + n) in TV 



Ari/3 



oo. 



Note that the factor {tn/tm)^'^~^^^'^ , which is irrelevant in calculating determinants, was omitted in 
the second line in the above equations. Put u = —A to obtain the expression (I5.13p . Similar evaluation 
in — > oo of (j4.25p can be done also for m > n. | 

The infinite system obtained by the soft-edge scaling limit (15. 7p is temporally homogeneous, but 
spatially inhomogeneous as shown by the "Airy kernel" (j5.13p . Prahofer and Spohn [78] and Johansson 
|46j studied the right-most path in the present system and called it the Airy process A{t). For a given 
t > 0, A(t) has distribution of the celebrated Tracy- Widom distribution, which is governed by the 
Painleve II equation [94J . Recently Tracy and Widom derived a system of partial differential equations 
(PDE), which govern the Airy kernel (j5.13p [95j . They also discussed other determinantal processes 
by PDE [96]. See also [DE]. 



6 DETERMINANTAL PROCESSES ASSOCIATED 
WITH SPECTRAL PROJECTIONS 

6.1 Spectral Projections 

First we note that, following the notations in Section 2, if we set r = logt,r' = logt',C = y/\/2i, C' = 
x/\/2F, and 'KN{t,y;s,x)dy = Kn{t,C',t' ,(')dC with d( = dy/\/2i, the matrix-kernel of the determi- 
nantal process of noncolliding BM with finite number of particles N < oo given by (I4.25P is rewritten 
as 

g K|e'''-''"'P^IC'> if t<t' 
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where Hip is given by ()2.12p . and V^p is a projection operator defined by 

V^= Yl (6-2) 



0<fc<Af-l 



It is called the extended Hermite kernel in 
The trigonometric functions of the form 



S-(VXx) = —=^-—sin(\/Xx), 



can be regarded as the generalized eigenfunctions of the Hamiltonian 

with spectrum A > 0. Here S_ is an odd function (parity p = — ) and 5+ is an even function (parity 
p = +), respectively. Consider an operator TYsin such that (x|?Ysin|y) = ^{x — y)Tlshi with (j6.3p and 
introduce a set of its eigenvectors {|A, p; sin) : A > 0, p = ±}; 

Wsin|A, p;sin) = A|A, p;sin). 

Since (x|A, p; sin) = (A; p; sin|a;) = 5p(\/Ax), we can confirm the completeness of the set 

/ (iA(x|A, p;sin)(A, p;sin|y) = — / du cos{u{x - y)) = 5{x - y). 
p^Jo ^'^ J~oo 

If we change the variable in the Airy differential equation (|5.9p by x ^ u — A, we have 

(f 

— r^Aiix - A) + xMix - A) = AAi(x - A). 

That is, the Hamiltonian 

^Ai = -^ + a: (6.4) 

has R as spectrum and the Airy functions of the form Ai(x — A) are its generalized eigenfunctions. 
We can consider the corresponding operator TYai and its eigenvectors {|A; Ai) : A G R}, 

WAi|A;Ai) = A|A;Ai), 

where (x|A; Ai) = (A; Ai|x) = Ai(x — A). We find the completeness 

/>oo 

dA(x|A;Ai)(A;Ai|y) = / d\ k\{x - X)Ki{y - \) = 5{x - y). 



Remark 3. The 2{v + l)-dimensional squared Bessel process (BESQ^) y(^)(t) is defined as a 
solution of the SDE 

dY^^^t) = 2\/Y(''){t)dB{t) + 2{iy + l)dt, u>-l, (6.5) 
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where B{t) is a one-dimensional standard BM |80lll8j . Its forward and backward Kolmogorov (Fokker- 
Planck) equations are given as 



d d'^ d 



(6.6) 



where — for the forward and + for the backward equations, respectively. We will see that the Laguerre 
polynomials {L'^{x)} play for BESQ,/ the similar role to the Hermite polynomials {Hk{x)} for the BM 
shown in Section 2. By considering the noncolliding system of BESQ,yS instead of BMs [59j, we can 
derive a finite determinantal process whose matrix-kernel is described using the orthonormal Laguerre 
functions (extended Laguerre kernel If we take the so-called hard-edge scaling limit, a spatially 

inhomogeneous infinite particle system is obtained, which is the determinantal process associated with 
the matrix-kernel 



dX e("-*)^ Ju (2^) 



J,{2^)^Jl{2^) - ^Jl{2^)J,{2^) 



1: 



y-x 

dAe-(*-")^J^(2v^)J^(2yA^) 



if t < s 



if t 



if t> s, 



(6.7) 



where J^(x) is the Bessel function, Jy{x) = Y^'^^^{-lY{x/2f''+''/{T{n + l)T{n+l + Ty)}, and J^( 
dJu{x)/dx. This kernel was called the extended Bessel kernel in [96]. See also 
see that J,y(2\/Ax) is the generalized eigenfunction of the Hamiltonian 



X 



. We can 



-Hj 



dx 5x ~'~ 4x' 



(6.^ 



with spectrum A > 0. We introduce the corresponding operator TLj and its eigenvectors {|A; J) : 
A > 0}, nj\X;J) = A|A; J), where {x\nj\y) = 6{x - y)Hj, (x|A; J) = (A; J|x) = Jy{2^/kic) with the 
completeness 

d\{x\\;J){X;J\y) = / dX M2VX^)M2y^) = 6{x - y). 
Jo 



Here we note the fact that the sine kernel K given by (jS.Sp . the Airy kernel IC given by ()5.13p . and 
the Bessel kernel given by (|6.7p are expressed in the same way as ()6.ip for the Hermite kernel IK at, 



K{t,y;s,x) 



. -{y\e 



if t<s 



l-V)\x) if t>s, 



if we assign the Hamiltonian 5{x — y)7i = {x\7i\y) and projection operator as follows instead of TC^ 
and (USD; 



(i) for the sine kernel, set 7i = "Hgin with (16. 3p and "Psin = 2. / 1-^' P' sin) (A, p; sin| 

.0 

(ii) for the Airy kernel, set H = Hm with dO]) and Vm = dX |A; Ai) (A; Ai|. 

J ~oo 

(iii) for the Bessel kernel, set 7i = 7ij with (16. 8p and Vj = dX\X; J) (A; J|. 

Jo 
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6.2 Effective Hamiltonians and Matrix-Kernels 



In the previous subsection, we claimed that the structure of the matrix-kernel of determinantal cor- 
relation functions is common both in finite particle systems and infinite particle systems. It should 
be noted, however, that even from the same finite system (e.g. noncolliding BM governed by the 
Hamiltonian Ti.^), depending on scaling limits, different kinds of infinite determinantal systems are 
derived, in which the Hamiltonian is replaced by appropriate effective Hamiltonians {e.g. TCsin and 
Ti-Ai) and spectral projection operator is modified. 

In the present subsection, we give a possible general consideration on the common structure of 
determinantal processes. First we note the following fact for a general form of effective Hamiltonian 
Tl = —a{x)d'^ /dx^ — h{x)d/dx — c(x), where a{x) ,h{x) , c{x) are sufficiently smooth functions with 
a{x) 7^ in an interval A C R. If we change the variable x ^ z following 

z{x)= r a{yr^/^dy, 







n^n = -d^/dz^-b{z)d/dz-ciz) with6(z(x)) = a{x)z"{x)+bix)z'{x) = a{x)-^/^{-a'{x)/2+b{x)}. 
Then we define 

r(z)=exp| — -y b{u)du^^ 

and if we perform a similarity transformation TL ^ TL = r^^TLr, the term of first derivative can be 
eliminated and we have the form TL = —d'^/dz'^ + liz). The transformation TL ^ TL is called the 
Liouville transformation. Then, without loss of generality, we can assume effective Hamiltonians of 
the form (the Sturm-Liouville operator ^93j ) 

^ = -^+9(^) (6-9) 

defined on A C R. 

Example 1. The effective Hamiltonians TLy,,TLsin,TLAi, and TLj are transformed to the form (|6.9|) 
with the following q{x), respectively 



—A':'-*). 0, I, [v^--]—, (6.10) 



where A = R for the first three cases and A = R+ = {x G R : a; > 0} for the last case. 

Let TL be the operator corresponding to (j6.9p : {x\TL\y) = S{x — y)TL, and define 

6t{x,y) = 6t{y,x) = (y|e-*«|x). (6.11) 
By definition, it solves the equation 

^^6tix,y) = -n6t{x,y) (6.12) 



with 



lim 6t{x, y) = 6 {x - y). (6.13) 
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Assume that H has a distribution uj{d\) of spectrum a{7i) = {A : 'H\\) = A|A)} with the complete set 
of eigenvectors {|A) : A G cr(7Y)}. Then the spectral representation of 6t{x,y) is given by 

5tix,y) = (y|e-*^ / ^ u;(dA)|A)(A|x) 

a{H) 

with ^x{x) = {x\X) = (A|x). We assume that 

^C>0 s.t. I dy{x -y)'^6t{x,y) < Cf, t e {0,1], e A, (6.14) 
Ja 

where C does not depend of t. 

Example 2. For the four examples ()6.10p . we have the following explicit expressions of 5t{x,y); 
q{x) = ^(x2-4) : 

^ oo 

5,{x,y) = ^E'/'"(f)'/'n(|)e-V2 

n=0 

= = exp < (x — v) coth - — -xy tanh - 

2v/vr(l-e-*) I 8^ 4 ^ V4 

X, y G A = R, 



/(x) = : 



/•oo ^ 

6t{x,y) = V / dA(y|e-*^-"|A,p;sin)(A,p;sin| 



fOO 

|x 

p=± 

— J --^e~'^*| cos(\/Ax) cos(\/Ay) + sin(\/Ax) sin 
1 2 

— / * cos(n(x — y)) 
71" Jo 

:exp<^- ^ ^/^ L x,yGA = R, 



/4^ I 4t 



(x) = X : 



/oo 
Ai(x - A)Ai(y - X)e'^^dX 
-oo 



1 f (x-y)2 t(x + y) , . „ 

:exp<^- ^ + — x,y£A = R, 



Airt 1 4i 2 12 



1 /""^ 

6tix,y) = - ^/xJu{\^x)^J^{VXy)e'^^dX 
^ Jo 



where Iu{z) is the modified Bessel function given by Iu{z) = "^^=0(2 / 2)'^'^^'^ / {n\T (v + n + 1)}. We 
can confirm that (I6.14p is satisfied in these four cases. 
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Now we consider a subset of spectrum cr{Ti.), (T-iJi) = {A G o'('H) : A < A*} with a specified level 
A* G cr(TC), and define the projection operator onto cr-iTC) 



so that 



and 



V = I ^ cj(dA)|A)(A|. 



Gt{x,y) = Gt{y,x) = {y\e'^V\x), 
gt{x,y) = g,{y,x) = -{y\e-''^{l-r)\x) 



K{x,y) = limgt{x,y) = limg^t{x,y) = {y\V\x), 



By definition 



p{x) = K (x , x) = {x\V\x) . 



Qtix, y) = g-tix, y) - 5t{x, y), 



d_ 

di 
d_ 

dt 



Gtix,y) = ngtix,y), 
Gt{x,y) = -ngt{x,y), 



and 



Moreover, by the completeness of : x G A}, Jj^dx\x){x\ = 1, we have the relations 



dy 5t{x,y)gt{y,x) = p{x), 



dy 5t{x, y)gt{y, z) = K{x, z). 



We define matrix-kernel K by 

K{t,y]s,x) 



gs-t{y,x) if t<s 
gt-s{y,x) if t>s, 



(6.15) 



(6.16) 

(6.17) 
(6.18) 

(6.19) 



(6.20) 
(6.21) 

(6.22) 
(6.23) 



(6.24) 



and consider the determinantal process H(t) = {JC(t)}, whose multitime correlation function is given 
by 

l<j<Nm,l<k<Nnfl<m,n<M 



(6.25) 



for any M > 0, any sequence {Nm}m=o °^ positive integers, and any series of observation times 
< ti < • • • < tM- 

The invariant measure of the process is the determinantal point field /i associated with K{x,y) 
given by (j6.17p . 
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The two-time correlation function p(0, {Xfn}',t, {?/„}) of the system is given by 
p{0,{xm};t,{y^}) = detM{t,yJxm), t > 0, m,n > 0, 



where 



M{t,yjxr, 



( p{xl) 


■■ K{xi,Xm) 


Gt{xi,yi) • 


• Qt{xi,yn)\ 


K{X2,XI,) ■ 


■■ K{x2,Xm) 


Gt{x2,yi) • 


■ Gt{x2,yn) 


K{Xm,Xi) ■ 


p{Xm) 


Qt{xm,yi) ■ 


■ Gt{Xm,yn) 


Qt{yi,xi) ■ 


■■ Qt{yi,Xm) 


p{yi) 


■ K{yi,yn) 




•• Gt{y2,Xm) 


K{y2,yi) • 


• K{y2,yn) 


\Gt{yn,Xi) ■ 


■■ Gt{yn,Xm) 


K{yn,yi) • 


PiVn) / 



(6.26) 



For a matrix A = {aij)i£jj£j with index sets /, J, we denote its submatrix as Ajiji = (aij)i^i'j^j' for 
I' C /, J' C J, and the complementary submatrix as = {0'ij)iei\i' ■ By using the relation 

(I6.19p . /3(0, {aim}; {j/n}) is expanded as 

p(0,{ic^};t,{y„,}) = detM(t 

mAn 

+ ^ ^ ^ (_l)^+ELi('ii+'»+60 

i=l l<ai< --<ae<m,l<bi<- -<bi<n 

X det D(t, yn\xm){a,}{b,} det M(t, 2/„|a;„)^'"+^^>^"^>, 



(6.27) 



where 



M{t,yjxr^ 



1 P{xi) 


■ K{xi,Xm) 


Gt{xi,yi) • 


• Gt{xi,yn)\ 


K{X2,XI,) ■ 


■ K{x2,Xm) 


Gt{x2,yi) • 


■ Gt{x2,yn) 


K{Xm,Xi) ■ 


p{Xm) 


Gtixm,yi) ■ 


■ Gt{Xm,yn) 


Q-t{yi,xi) • 


■ G-t{yi,Xm) 


p{yi) 


■ K{yi,yn) 


G~t{y2,xi) • 


■ G^t{y2,Xm) 


K{y2,yi) ■ 


■ K{y2,yn) 


\G-t{yn,Xi) ■ 


■ G-t{yn,Xm) 


K{yn,yi) ■ 


p{yn) / 



J:>{t,y„\xr. 



St{xi,yj)^ 



l<i<m,l<j <n 



It was shown by Shirai and Takahashi |84j that the Palm measure //^ coincides with the determi- 
nantal point field associated with the kernel defined by 



K'ix,y) 



1 



K{x,y) K{x,z) 
Kiz,z) ''^\K{z,y) K{z,z) 



det 



Note that 



For /,<7GCo(R), let 



P'{{Xm}) 



P{z) 



pi{Xm}U{z}). 



if, = ^ fix) and {g,i]) =J29{y), 



(6.28) 
(6.29) 
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and set a, (3 £ C. Define -F(^) = e"^-^''>\ G(^) = e^^^'^K Then the two-time generating function is 
defined by 

(6.30) 



E 



for t > 0. Let 
then we can see 



<^t(/,5;a,/3) = E[F(S(0))G(H(t)) 



oo oo ^ „ „ m n 

7 CHARACTERIZATION OF DETERMINANTAL PROCESSES 

In the previous section, we introduced a class of determinantal processes associated with spectral 
projections defined by effective Hamiltonians. Here we give properties of determinantal processes of 
this class, which can be derived from the common structure of correlation functions. 

7.1 Continuity 

Let C^(R) be the set of all infinitely differentiable real functions with compact supports and set 
{/,(,) ^ Ylxe^fi^) f ^ C'q°(R),^ G X. By a criterion of Kolmogorov (see, for example. [171 [12] ). 
the following lemma implies that X]xgS{t) /(^) continuous in time with probability one for any 
/ G Co°(R). Since X is separable with the vague topology, we can choose a countable set {fj}'^i C 
C^(R) such that — *■ C iii ~^ on j£ if and only if fj) fj) in n — > oo, Vj > 1. Then it 

implies that the determinantal process H(t) is continuous in the vague topology with probability one. 
That is, if the condition (|6.14p is satisfied, it can be expressed by 



it) 



with some real- valued continuous processes Xj{t),j £ N. 

Lemma 7.1 Let S(t) be the determinantal process, whose multitime correlation functions are given 
by 116. 25\) with 116.16]) and \6.24^ associated with an effective Hamiltonian \6. y\) defined on A C R. 
Assume that {6.14-^ is satisfied. Then for any f £ C^(R) 



E 



(/,H(t))-(/,H(0)) 



< c^^ tG(o,i], 



(7.1) 



where C does not depend on t. 



Remark 4. Theorems 15.11 and 15.21 are the limit theorems of the processes 
H^"''^(t) EE {Xi(iV + 2t), ■ ■ ■,Xn{N + 2t)} =^ E^'^it) 

E''^^''{t)^[Xi{N'/' + 2t)-aN{t),---,Xr,{N'/' + 2t)-aN{t)} ^ H^XO, 

with a^it) defined by (|5.11|) . in the sense of finite-dimensional distributions, where the multitime 
correlation functions of H'^™(t) and E^\t) are given by psin of ()5.4p and pAi of (I5.12p . respectively. When 
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we apply this lemma to the A^-particle system of noncolliding BM, which is a determinantal 

process associated with the Hermite kernel K^r given by (j6.ip . C depends on N. Careful estimation 
gives an upper bound on C, which is uniform in N, and tightness is established both in {H^"^(t)}ArgN 

and {H^^°(t)}ArgN- 

Proof of Lemma 7. 1 . Here we use the notation 

Pm,n{p^m+n) — {^1 ; ' ' ' ; i ; • • • ; ^m+n})' 

The left-hand-side of (17. ip is given by d^^t{—f, f] a)/d^a\a=o from the two-time generating function 
given by (|6.3Up and (j6.3ip . which equals 

r ^ 

/ dx4 /(a;i)|p4,o(a^4) - 4^3,1(0:4) + 6^2,2(3^4) - 4^1,3(2:4) + po,4(a;4)| 
■^^^ i=i 

» 3 3 

+ / dxsl\f{xi){2y2f{xj)ps,o{xs) + 6{-f{xi)-f{x2) + f{xs))p2,i{x3) 

r=i ^ U 

3 

+6(/(xi) - f{x2) - f{xs))pi,2ix,) + 2 ^ f{xj)po,3{x,)] 

j=l 



+ 



dX2 n f{xi)[{2f{xif + 3/(xi)/(x2) + 2f{x2f){p2fl{x2) + P0,2{X2)) 
i=l 

(-4/(xi)2 + Qf{x^)f{x2) - ^f[x2f)pi,l{x2)} 



+ 



+ I dxi f{xi)'^^^pifi{xi) + po,i{xi)y 



Since P3,o{xi, X2, X3) = po,3ixi,X2,X3), p2,i{xi,X2,X3) = pi,2{x3,X2,xi), etc, the above quantity is 
twice of 



j ^ n/(a^i){p4,o(a^4) - 4/93,1(334) + 3/92,2(a;4)| 

3 3 

dX3 n fi^^){^Yl fi^j)pM^3) + 6(-/(xi) - f{x2) + f{x3))p2,l{^3)} 
i=l j=l 

+ / dX2T\f{x,)\i2fixif + 3f{xi)fix2)+2fix2f)p2,oix2) 

+(-2/(xi)2 + 3/(xi)/(x2) - 2/(x2)')pi,i(a;2)} 
+ / dxi /(xi)^/3i,o(xi). 



We put Mm,n{xm+n) = M{t,Xm+i, ■ ■ ■ ,Xm+n\xi, . . .,Xm), D(a;„) = D(t,a;„|a;„) and 

Di, 1(332) = -St{x2,xi)gt{xi,X2), 



D2,i(3=3) = ^(-l)^+=^5t(x3,x,)detM2,i(a;3)^=^>« 
1=1 
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D3,i(a;4) = J](-l)*+'5i(x4,x,)detM3,i(a;4)^^>«, 



i=l 
2 4 



Prom (|6.27p we have 



D2,2(a;4) =Y,Y.(-'^y^'^t{xj,Xi)detM2,2{xA)^'^^'\ 

i=l j=3 

02,2(2:4) =detD(a;4){3,4}{i,2} det M2,2(a;4)^^'^"^'^^ 



Pifii^i) = detMi,o(a;i), 

Pi, 1(3^2) = detMi,i(iC2) - Di,i(a;2), 

^2,1(333) = detM2,i(a;3) - 02,1(0:3), 

P3,i(a^4) = detM3,i(a;4) - 03,1(034), 

P2,2{^4) = detM2,2(a;4) -02,2(0^4) +D2,2(a;4)- 



We divide / into four terms / = X^j=i Ij with 

h ^ [ dxi /(xi)^detMi,o(xi) 
Jk 

2 

dx2 n f{xi){-2f{xif + 3/(xi)/(x2) - 2/(x2)')Di,i(o;2), 

i=l 

2 

dx2 n f{xi)[{2f{xif + 3/(xi)/(x2) + 2/(x2)=')det M2,o(o;2) 

i=l 

+{-2f{xif + 3/(xi)/(x2) - 2/(x2)')detMi,i(o;2)} 
3 

- / dx:, J]/(xi)6(-/(xi)-/(x2) + /(x3))D2,i(o;3) 
i=i 

r 4 
+ / c^a;4 JJ/(xi)3D2,2(a^4), 

i=i 

3 3 

dx^ n /(^0{2 /(a;,)det M3,o(o;3) + 6(-/(xi) - /(X2) + /(x3))det M2, 
i=i j=i 

- / dx^ n^(^i){ -^^3,1 (3^4) +302,2(0:4)}, 
4 

'2^a;4 JJ/(xi)|detM4,o(a;4) - 4det M3,i(o;4) + 3det M2,2(o;4)}- 

1=1 

Using (|6.22p we have 

h = -2 / f{xi)f{x2fGt{xi,X2)5t{xi,X2)dxidx2 

JA2 



-2 / /(Xi) f{x2)Gt{xi,X2)dt{xi,X2)dxidX2 

Ia2 
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A2 



f{xi)'^f{x2fGt{xi,X2)6t{xi,X2)dxidx2+ / /(xi)^/)(xi)dxi 



1 ^ 

f{xi) - f{x2)j Gt{xi,X2)6t{xi,X2)dXidX2 



By the assumption (I6.14p we have 



h<Cit\ tG(0,l], 



(7.2) 



where Ci does not depend on t. 

Since Gt{x,y) = K{x,y) + dQt{x,y)/dt\t=o t + 0{t^), we have for any k,£ 



Then 
where 



detMk,i{xk+e) = det Mfe_i,£+i(a;fc+£) + 0{t'^). 
Ij=ij + 0{t^), J = 2, 3, 4, 



r ^ 

/ dx2 Ylf{xi)6f{xi)f{x2)detM2.o{x2) 

3 4 

- / dX3 n/(xi)6(-/(xi)-/(x2) + /(x3))D2,l(cC3)+ / dx^ llfiXi)3B2,2{Xi), 

r ^ 

- / da;4 T\f{xi){-4D3,i{xi)+3B2,2ixi)}, 



and /4 = 0. Since the estimate ()7.2p was obtained, it is enough to show that 

Ij<Cjt^ iG(0,l], j = 2,3 
for the proof, where Cj's do not depend on t. Using eqs. (I6.22|) and (I6.23|) . we obtain 



(7.3) 



/ dX4 

JA4 



Gt{xA,Xl) Qt{Xi,X2) 
Gt{x3,Xl) Gt{x3,X2) 



^fixiffix2f6tixi,X4)5tix2,X3) 



+3/(xi)/(x2)/(x3)/(x4) 
+6/(xi)/(x2)/(x3)(-/(xi) - /(X2) + /(X3)) 

dx4 



A4 



Gt{xA,Xi) Gt{x4,X2) 
Gt{x3,Xi) Gt{x3,X2) 



6t{x4:,Xi) St{x4,X2) 
(5t(x3,Xi) St{x3,X2) 

(5t(x4,Xi) 5t{x4,X2) 

(5i(x3,Xi) 5t{x3,X2) 

St{x4,Xi) (5t(x4,X2) 
St{x3,Xi) dt{x3,X2) 



3Fix^), 



where 



F{X4) = /(Xi)/(X2){/(X1)/(X2) + /(X3)/(X4) " 2/(xi)/(x3) - 2/(x2)/(x3) + 2f{x3f]. 

By simple calculation we see that 

F(xi, X2, X3, X4) + F(xi, X2, X4, X3) + F(x3, X4, Xi, X2) + F(x3, X4, X2, Xi) 
= -2(/(xi) - /(X3))(/(X3) - /(X2))(/(X2) - /(X4))(/(X4) " /(xi)). 
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Then 



dxA 



A4 



Gt{x3,Xi) gt{x3,X2) 



Stix4^,Xi) 6tix4,X2) 
(5i(x3,Xi) 6t{x3,X2) 



A4 



dX4 5tix4, Xi)5tix3, X2)ifixi) - /(X4))(/(X2) - /(xa)) 



gtix4,Xl) Qt{x4,X2) 
Qt{x3,Xi) Qt{x3,X2) 



A4 



(/(Xl)-/(X3))(/(X2)-/(X4)) 
dXi 5t{x4,Xi)5t{x3,X2){f{xi) - f{Xi)){f{x2) - fixs)) 



x[{Gt{x4,Xl)gt{x3,X2) - gt{xi,X2)gt{x3,X4)){f{xi) - f{xi)){f{x2) - f {xz)) 
-Qt{xA,Xi)Qt{x3,X2){f{xi)f{Xi) + /(X2)/(X3)) 
+Qt{x2,Xi)Qt{xi,X3){f{xi)f{x4) + /(X2)/(X3))| 

/ dxi 5t{xi,xi)5t{x3,X2){f{xi) - f{xi)){f{x2) - fixs)) 

x\^iGtix4,Xi)gtix3,X2) - gtixi,X2)gtix3,X4))ifixi) - f{xi)){f{x2) - f {xz)) 

+{Qt{x2,X4)Qt{xi,X3) - Qt{x2,xi)gtix4,X3)){f{xi)f{x4,) + /(a;2)/(a;3) ) | 

dX4 6t{x4,Xi)6t{x3,X2){f{xi) - f(x4))(f{x2) - /(xs)) 



A4 



{0t{x4„Xi)Qt{x3,X2) - Qt{xi,X2)Gt{x3,Xi)){f{xi) - fixi)){f{x2) - /(X3)) 
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+ ( ^^gt{x^,X2)^^gt{x^,X2) - 0t{x2,x,) ^^^^^^ 



Gt{xi,X2) 



X(/(X1)/(X4) + /(X2)/(X3))(X4 " Xi)(x3 " Xs)] + O(t'). 

By the assumption ()6.14p we obtain (|7.3p for j = 2. 
Since 

(-l)^+^5t(xi, X4) det M3,l(xi, X2, X3, X4)^^>^^> = (" 1)^+^5* (X2 , X4) det M3,1(X2, XI, X3, X4)^'>^'^ 

= {-lf+^6t{x3,X^) det M3,1(X3, X2, XI, X4)^3^^^>, 
(-l)^+^5t(xi, X4) det M2,2(X1, X2, X3, X4)^^>^^^ = (- l)^+=^5t (xi , X3) det M2,2(X1, X2, X4, X3)^^>^'^ 

= {-lf+Ht{x2, X3) det M2,2(X2, XI, X4, X3)^2^^^^ 

= (-1)2+4 5i(x2, X4) M2,2(X2, XI, X3, X4)^'^^^>, 



we have 



12 / dxi T\ f{xi)5t{xi,X4) 



K{x2,Xi) p{x2) K{X2,X3) 
K{x3,Xi) K{X3,X2) P{X3) 
Qt{Xi,Xi) Qt{xA,X2) Qt{xA,X3) 



+ 



K(x2,Xi) p(x2) ^_t(x2,X3) 
^i(x3,Xi) ^t(x3,X2) P(X3) 
Qt{xA,Xi) Qt{Xi,X2) K{xi,X3) 



Then we have ()7.3p for j = 3. This completes the proof. 
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7.2 Bilinear Forms 



Since the Fredholm determinant of the dual operator coincides with that of the original operator, the 
reversibility (and also the stationarity) of the present determinantal processes is guaranteed. 

Osada constructed X-valued reversible processes, which have determinantal point fields as their 
reversible stationary measures by Dirichlet form approach. With the Palm measure /i^ the Dirichlet 
form of Osada can be written as 



£{F, G) = dzp{z) |^^^(dr?)^F({z} U ^)^G({z} U r?) 



(7.4) 



for local smooth functions F, G [76] (see also [88] ) . By the general theory of Dirichlet forms [38] , his 
processes are diffusion processes (i.e., continuous strong Markov processes). For our class of deter- 
minantal processes, we found the following fact (Proposition 17. 2( ). It suggests that our determinantal 
processes are identified with Osada's processes. 

A function F on the configuration space X is said to be polynomial, if it is written of the 
form F{Q = F{{fi, OAf 2,0^" ' Afk^O) with a polynomial function F on Il^,k G Nq, where 
fj G C^(R),1 < j < k. Let p be the set of all polynomial functions on X, which is a dense 
subset of L^(jC, /i) ; the space of square integrable functions on X with the determinantal point field 

Proposition 7.2 Let H(t) be the determinantal process, whose multitime correlation functions are 
given by i6.25\) . Then for F,G £ p we have 



dt 

where E is given by |7.^|j. 



|^E[F(H(0))G(H(t))]^^^ = ^:(F,G), (7.5) 



Markov property of determinantal processes has been studied by Borodin and Olshanski [14^ I16j. But 
we can not prove that the infinite determinantal processes in our class are Markovian. The equivalence 
between our processes and Osada's is not yet established. 

To show this proposition, it is enough to consider the case that F and G are of the form 

F(^) = exp I ^ OLm'^fmix) I , G{i) = exp I ^ /3m ^5m(a 

ym=l I'Gg I \rn=\ xe£, 

with k £ No, am, (3m G R, /m,5m G C^{R), l<m<k. Then if we set xii^) = exp {Ylm=i (^rnfn 
1 and X2{x) = exp ( Em=l Pm9mix)) - 1, 



m = 11(1 + XI (x)), GiO = Uil + X2{x)). 



The left-hand-side of ()7.5p equals 



-g^Mifm}, {gm}; {am}, {Pm})\^^^ 



/ dy Y[xi{xj)llx2{yk)Trp(.0,{xm};t,{yJ) 

m=0n=0"'-'^- 7=1 k=l 



t=0 
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Since the Palm measure is a determinantal point field associated with p^{{xm}) (see (I6.28P ) the 
right-hand-side of (|7.5|) equals 

= I dzp{z) ^'^^^^^ ^^^^^^ V [ dXm ']\{xi{xj)X2{Xj)+Xl{Xj)+X2{xj)\p''{{x.ra]). 

oz oz m. Ja™ -^^ ^ ' 

Hence Proposition 17.21 can be derived from the following lemma. 
Lemma 7.3 For any 171,71 £ IS we have 



lim dXm dy^ Y\xiixj)Y[x2{yk)w:PiO,{xm};t,{y„}) 

t^OJA'^ J An dt 



\^ mini f dxi{z)dx2{z) 

^ (m-^)!(n- £)!(£- l)!y/"^^"^ 



dz dz 



„ /•I- m—£ n~£ £—1 

X dXm-£ dy^__^ dwi_i T[ Xi{xj)T[ X2{yk)T\ Xi{wi)X2{wi 

J^-' 7a"-^ 7a^-i f^^ ^J-^ t\ 

xp^{{xm-e} U U {wi_i}). 



Proof of Lemma 7.3. We use the expansion formula (j6.27p of the two-time correlation function to 



calculate the integral 



» „ m n 

1=1 dXm / dy^Y[xi{xi)T\x2{yj)piO,{x.,n};t,{yn}). 
By permutation invariance of the integrand we put ai = m — £ + i,bi = i with i = 1,2, . . . and then 



dXm dy„ M xi(xi) M X2(yi)detM(i 
""^^ fm\ (n\ r r " 

+ dx^ dy^ n u ^^(y.) 

X det D(t, yn\xm){a,}{b,} det M{t, yJ^„,)^'"+^^>^"^> 

„ m n 

dxrn I dyn\\xi{xi)\\x2{yj)<^eim{t 

•^A i=l J=l 

where Mf(a;™;y„) = M(t, 2/„|a;„){"^+i'-'"^+^H™-^+i,-,™}. By using (f6J3]) . dHZ]), ((6:20]) and (lOTD . 
we see that 



I A" 

mAn 



l^'^lu dXm dynT\xi{xi)Y\x2{yj)detM{t,yJxm) = 0. 

i^U Or / Am /An 



(7.6) 
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From (f6J2]) and IK2Qh 

•^m j Un) 

■/A™ -/A" j,^-^ 

= X! / '^^"i / c^y^n^i^'^^) n^2(2;j)^<5f(a;™_£+p,?/p) J| (5t(a;m_f+fe, y^) det (a;™; y„) 

p^i^A'" Ja" fe=i,fc5^p 

m—i ^ ^ m n £ 

— 1 J A™ J A" K_i 



p=l 



X det 



K{Xm,Xl) 

G-t{ye+i,xi) 
Q-t{yt+2,xi) 

\ Q -t{y n,Xl) 



i=l j=l k=l 






dyt + 2^-t{ye+2,Xp) 

9LG-t{yn,Xp) 



K{xi,Xm-i) Gtixi,yi) 

K{x2,Xm-i) Gt{x2,yi) 

K{x ) Gt{xm,yi) 

G-t{yi+i,Xm-e) K{yi+i,yi) 

G~t{ye+2,xm-e) K{yi+2,yi) 

G-t{yn, X m—i ) K{yn,yi) 



n „ „ m n £ 

^/ dxm dy„Y[xiixt)'[[x2{yj)'[[Stixm-e+k,yk) 

1 -/a™ ja" ._i , 1 



p=l 



X det 



K{xi,xi) 

K{x2,Xi) 

K{x„i,xi) 
G-t{ye+i,xi) 
G~t{ye+2,xi) 



i=l j = l k=l 

K{xi,Xm~t) Gt{xi,yi) 

K{x2,X„i-£) Gt{x2,yi) 
K{x,n,X„i-i) Gt{Xm.,yi) 

G-t{ye+i,Xm~i) K{ye+i,yi) 
G-t{ye+2,Xm-i) K{ye+2,yi) 



\ G-t{yn,Xl) 

hit) + hit) -hit), 



dlGt{xi,yp) 



dyGt{x2,yp) 



dy^GtiXm,yp) 








Gt{xi,yn) \ 

Gt{x2,yn) 
Gt{Xm,yn) 

K{yi+i,yn) 
K(yt+2,ya) 

K{yn,yn) ) 



Gt{xi,yn) \ 

Gt{x2,yn) 
Gt{Xm,yn) 

Kiyi+i,yn) 
K{yi+2,yn) 



K{yn,y,i) J 



G-t{yn,Xm-i) K{yn,yi) 

where we have used the abbreviation dy = d"^ /dy^. By partial integration 
h{t) = X! / / II X2{yj)\\_6t{x„i-i+k,yk)^\x2{yp)<^G^M.l{x^;y^)\ 

Ja- JA" i=l J5^p ^ ^ 

f. „ „ m n I q2 

= X! / / dy^'^xi{xi)'Wx2{yj)'W5t{x,n-i+k,yk)-^'^QtM.l{x^]y^) 

p=i-J^"' -^A" oyp 

£ „ m n £ 

p^l JA™ JA" . , . , 



(7.7) 



dyl 

hl{t)+h2{t). 



det M*(a;„;y„) + 2 



9x2(yp)9detM^( 



9yp 



5yp 



By simple calculation 
hi{t) + hit) - h{t) 

^/ da;„ / dy^'[Yxi{xi)Wx2{y])W5t{xm~i>+k,yk) 

1 •/A"» JA" ._i „_i , 1 



p=l- 



fc=l 
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X det 



p=i 



X det 



1 A(Xl,Xl) 


K{Xi,Xm-e.) 


Gt{xi,yi) 




U 


Gt\x\,yn) \ 




K{X2, Xi) 


K{X2, Xm-e) 


Gt{x2,yi) 







Gt{x2, y-n) 




K{Xm,Xi) 


Ki^Xfyi. Xjyi — ^^ 


Gtixm-.yi) ■■■ 







Gt{Xm,yn) 




G-t(ye+i,xi) 


■■ G-t{yi+l,Xm-l) 


K(yi+i,yi) ■■■ 


dy^K{yi+i,yp) ■■■ 


K{yi+i,yn) 




G-t{yi+2,xi) 


■■ G-t{yi+2,Xm-l) 


K(yi+2,yi) ■■■ 


dy^K{yi+2,yp) ■■■ 


K{yi+2,yn) 




\ G-t{yn,Xl) 


G-t{yn,X,n-e) 


K(yn,yi) 


dy^K{yn,yp) ■■■ 


K{yn,yn) J 




/ dXm / C?y„ 
J A™ JA" 


n^l(^*) n^2(yj) W St{Xrn-i+k,yk) 










/ A(Xi,Xi) 







-i) 


Gt{xi,yi) 


Gt{xi,yn) \ 


K{x2,Xi) 





K{x2,Xm 


-i) 


Gt{x2,yi) 


Gt{x2,yn) 


K{Xm,Xi) 





KijCjYi , Xui 


-i) 


Gt{xm,yi) 


Gt{Xm,yn 


) 


G t (vp-s-^ ■ Xi ) 


•• d^t+i^-tiye+i,xp 


) ■■■ G-t{yi+i,x 


m-l) 


K{ye+i,yi) 


■■■ K{yi+i,y 


n) 


G-t{ye+2,xi) 


■■ 9l^^^G-t{ye+2,Xp 


) ■■■ G-t{ye+2,x 


vn—i) 


K{yt+2,yi) 


■■■ K{yi+2,yn) 


\ G-t{yn,Xl) 


9y„G-t(yn,Xp) 


G-t{yn,Xm-f) 


K{yn,yi) 


K{yn,yn 


) ) 



^ / dx^ j dy^'^xi{xi)'Wx2{yj)'W5t{x,n-i+k,yk) 

3=^+1 i = l j = l k=l 

( K{xi,xi) ■■■ K{xi,Xm-t) Gt{xi,yi) 

K{x2,Xi) ■■■ K{x2,Xm-t) Gt{x2,yi) 



X det 



K{Xra,Xi) 

G-t(ye+i,xi) 
G-t(yi+2,xi) 

V G-t{yn,Xl) 



-^{Xjji^ Xyn — l) 

G-tiye+i,Xm-e) 



Gt{x„i,yi) 
K(yi+i,yi) 



G-t{ye+2,Xm-i) K{ye+2,yi) 
G-t{yn,x„,-i) K{yn,yi) 



dy^Gtixi,yp) 
dl^Gt{x2,yp) 

dl^Gt(xm,yp) 





Gt{xi,yn) \ 

Gt(x2,yn) 
Gt{Xm,yn) 

K{yi+i,yn) 
K{yi+2,yn) 

K{yn,yn) J 



By using ()6.13p and (I6.21|) we have 



lim<j/n(t) + /2(t)-/3(t)l>=0. 



(7.8) 



Suppose that Xm-i+i — *'yji^ = l>2,.... 



dyp 

for any p = 1,2, ... ,i. Hence we have 



det Mi{xm;yn) p{{xm-e} U {yj), 

ddetMiiXralVnl _^ 1 dp{{x,rn-£} U {Vn}) 

^ 2 



dyp 



in t ^ 



f „ „ m—i £ n 

hm/i2(t) = ^2 dxm^i dy^Ylxiixi) TT Xi{yk)X2{yk) TT X2{yj) 

p=l-^^^ -^^^ t=l k = l,k^p j=i+l 



dX2{yp) dp{{x,n-i] U {?/„}) 



dyp 



V/ dxm-i dy^Wxiixi) TT xi{yk)x2{yk) TT X2{yj) 

^ ' Km-t /An 

p=l-^'^ 1=1 k=l,k^p j=i+l 



37 



m—l 



n n III, t, lb 

V/ dxm-e dy^Tlxiixi) TT Xi{Vk)X2{yk) TT X2{yj 
p=l'''^ '''^ i=l fc=l,fc^p j=i+l 

^dxiiVp) dx2{yp) 



dyp dyp 



-pi{x^^l}\j{yj) 



m—l 



dxm~t dz dwi^i / dyn-eT\xi{xi)T\xi{yk)X2{yk)Y\x2iwj) 
pi{xm-i} u {z} u {wi_i} u {y„_ J) 



5x1 (2) 9x2(2; 



dz dz 

= / ^^^^^'^ 9X2iz) I dx,n~e I dwe_i [ dyn_^ 

Ja OZ dz Ja^-I JaI-I Jl^n-l 

m—£ £—1 n 

X n n xi{yk)x2{yk) n X2K)p^({a;„-4 u {w^-i} u (7.9) 
i=i fc=i j=i 

Combining ^L^ . ([TT]) . ([TIH]) and ([73]), we obtain Lemma El | 

8 CONCLUDING REMARKS 

The study on noncolliding BM and determinantal processes reported in this paper will be extended 
in several directions. We would like to give some of future problems below. 

(1) In Section 4 we let the initial configuration z/q be the GUE-eigenvalue distribution, and shown 
that the system is determinantal. If we let uq be the eigenvalue distribution in the Gaussian 
orthogonal ensemble (GOE), 



N 



(j3.1ip becomes 

ut{x) oc hN{x) / fN{t -to,x\x^^'>)dx^'^'^ 

= -^hNix) [ /^(t-to,a;|a;(°))sgn(/i;v(a=^°^))dx(°). 
Instead of the Heine identity (14. lip , we should use the de Bruin identity |25j 



•2) 



Wa 



dx det 

l<j,k<N L' 



9j{xk) 



Pf 



l<j,k<N 



R JR 



dx / dx sgn{x — x)gj{x)gk{x) 



for integrable continuous functions gj,l < j < N. As shown in Appendix A of [57J, for example, 
the generating function of multitime correlation functions is then expressed by the Fredholm 
Pfaffian [73] and the system becomes a Pfaffian process, in the sense that any multitime cor- 
relation function is given by a Pfaffian. Such Pfaffian processes have been studied by many 
authors [13 IMliSTl dZl [H21I31 Hi]- The systems studied in [33 EH [Ml E] are also Pfaffian 
processes, since the 'quaternion determinantal expressions' of correlation functions, introduced 
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and developed by Dyson, Mehta, Forrester, and Nagao [SU [651 EHl ED EH], are readily trans- 
formed to Pfaffian expressions. As implied by Proposition \'S.2\ the system exhibits a transition 
from GOE distribution to GUE distribution [671 [TTJ [52] . Continuity of sample paths and general 
characterization of infinite Pfaffian processes will be interesting problems. The case with other 
initial distribution (in particular, when it has continuous parameters) will be interesting j55ll57j. 

(2) As explained in Sections 1 and 3, the present noncolliding BM is the /i-transform of the 
absorbing BM in the Weyl chamber (jl.lOp of type Aj^^i. We can find appropriate /i-transforms 
of the absorbing BMs in the Weyl chambers of types Cat and Dat- The obtained noncolliding 
diffusion processes are stochastic versions of non-standard random matrix ensembles, which were 
called the class C and class D, respectively, by Altland and Zirnbauer [3l[31|99]. The stochastic 
version of the chiral GUE, realized by the noncolliding squared Bessel process, was studied 
by Konig and O'Connell [59], which is also obtained as an /i-transform of the absorbing BM 
in the Weyl chamber of type Cn- See [HI [55] for more details. Systematic classifications of 
determinantal and Pfaffian processes will be important. 

(3) There are many other examples of finite and infinite determinantal processes, which are not 
considered in the present paper. Markov processes on partitions (Young diagrams) have been 
studied and determinantal processes associated with other types of projections than ours have 
been reported [731 [HI [E]- The determinantal processes, whose kernels are expressed using 
multiple orthogonal functions {e.g. the Pearcey kernel), are discussed in [201 [HI [121 [23 151 
[Ml [971 [Ml [M]- The consideration given in Sections 6 and 7 in the present paper should be 
generalized. 
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